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ABSTRACT 



Recent observations of luminous Type Iln supernovae (SNe) provide com- 



pelling evidence that massive circumstellar shells surround their progenitors. 



' In this paper we investigate how the properties of such shells influence the SN 



lightcurve by conducting numerical simulations of the interaction between an 
expanding SN and a circumstellar shell ejected a few years prior to core col- 
lapse. Our parameter study explores how the emergent luminosity depends on 
a range of circumstellar shell masses, velocities, geometries, and wind mass- 
loss rates, as well as variations in the SN mass and energy. We find that the 
shell mass is the most important parameter, in the sense that higher shell 
masses (or higher ratios of M s /j e ;//Ms7v) lead to higher peak luminosities and 
higher efficiencies in converting shock energy into visual light. Lower mass 
shells can also cause high peak luminosities if the shell is slow or if the SN 
ejecta are very fast, but only for a short time. Sustaining a high luminosity 
for durations of more than 100 d requires massive circumstellar shells of order 
10 Mq or more. This reaffirms previous comparisons between pre-SN shells 
and shells produced by giant eruptions of luminous blue variables (LBVs), 
although the physical mechanism responsible for these outbursts remains un- 
certain. The lightcurve shape and observed shell velocity can help diagnose 
the approximate size and density of the circumstellar shell, and it may be pos- 
sible to distinguish between spherical and bipolar shells with multiwavelength 
lightcurves. These models are merely illustrative. One can, of course, achieve 
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even higher luminosities and longer duration light curves from interaction by 
increasing the explosion energy and shell mass beyond values adopted here. 

Key words: hydrodynamics — methods: numerical — stars: mass loss — 
stars: supernovae (general) — stars: winds, outflows 



1 INTRODUCTION 



The luminosity of a supernova (SN) results from energy in put by a com bina- 



tion of radioactive decay and shock kinetic energy (see e.g. 



Arnett 



1996), and 



for a Type II SN, the shape of the light curve depends on quantities l ike the 



star's initial radius, ejecta ma ss, and explosion energy (jArnett 



2004; 



Kasen fc Wooslev 



1996; 



Young 



20091 ) . For SNe with small initial radii, like SNe of 



Types la, lb, Ic, and peculiar SNe II like SN 1987A that result from blue su- 
pergiants, most of the shock-deposited thermal energy imparted to the stellar 
envelope is converted back into kinetic energy through adiabatic expansion, 
so nearly all of the observed luminosity comes from the radioactive decay of 
56 Ni and 56 Co. In "normal" SNe II-P that result from the explosions of red 
supergiants (RSGs), however, the large initial radius allows some modest frac- 
tion (typically 1-2%) of the shock-deposited thermal energy to be radiated 
away, powering much of the plateau of the lightcurve, although the vast ma- 
jority still goes into expansion energy. At late t imes, even SN e II-P have their 



Hamuv 



20031 ). 



luminosity powered by radioactive decay (e.g 
Subsequently, however, as the fast SN ejecta expand, they can collide with 
dense circumstcllar or interstellar material (CSM/ISM) that may surround 
the SN. As a result, additional kinetic energy may be transformed once again 
back into thermal energy through shock heating, which in turn may be lost by 
radia tive cooling if a dense radiative shock forms (e.g. 



Chevalier fc Fransson 



2008). This can enhance the luminosity for long after the explosion: Som 



SNe remain radio luminous for decades ( Montes et al 



2002 



Van Dvk et al 



1998; 



19931 ). and this interaction may power a visible 



nova remnant (SNR) su ch as Cas A for hundreds of years (|Chevalier 



; Williams et al. 


a visible s 


uper- 


( Chevalier 


1977; 



Chevalier fc Oishi 



2003). On the other hand, if the collision with dense CSM 



happens immediately after the explosion, it may significantly alter the spec- 
trum and light curve of the SN itself. This latter scenario is gen erally thought 



to be the case for the observed sub-class of Type Iln supernovae (jSchlegel 



Filippenko 



1990; 



19971 ). where the "n" corresponds to "narrow" or intermediatc- 



width H lines from the shock-heated C SM gas or decelerated SN ejecta (e.g. 



Chugai fc Danziger 



1994; 



Chugai 



2001). 
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In a normal SN, the expected results of radiative cooling and reheating of 
the SN ejecta due to radioactive decay yield can be estimated from analyt ical 
models of stellar structure and explosion physics ( Matzner fc McKeelll999 ). In 
SNe with strong CSM interaction such as the observed class of Type Iln SNe, 
however, the effects of collisions between an expanding SN and its circumstel- 
lar gas are harder to predict with ab initio calculations. They depend highly 
on the density and morphology of the CSM, which in turn depend on the 
unknown mass-loss behavior of the star in the few years prior to core collapse 
— potentially different for each object. A wide variety of CSM environments 
are possible, leading to a wide diversity of observed lightcurves and spectral 
properties. 

Recent observations of luminous T ype Iln supern ovae such as SN 2 006gy 



( Smith et al 



2007 



Ofek et al 



2008a) have 



20071) and SN 2006tf ([Smith et al. 
stretched the boundaries of our understanding of SNe Iln. Their extreme 
luminosities yield strong evidence that the progenitors of these SNe were sur- 
rounded by massive shells, presu mably ejected in precursor eruptions during 



the final years of stellar evolution ( 


Smith et al.ho07 


2008a. 


2010 


Smith & McCrav 


2007 


Wooslev. Blinnikov & HeserLoOTL 


Smith et al. 


(200 


7 ) pointed out that 



the physical properties (mass, speed, H composition) of these mass ejections 
were analogous to those observed for giant eruptions of luminous blue variables 
(LBVs), and espec ially reminiscent of the giant 1843 eruption of r\ Carinae 



(Smith et al 



2003). As the SN ejecta expand, they collide with the recently 
ejected CSM shell and this collision significantly decelerates the SN expansion, 
transforming kinetic energy back into thermal energy at the collision front, 
producing a brilliant fireworks display. The remarkably high luminosity and 
long duration of the observed emission from SNe 2006gy and 2006tf imply 
that the circumstellar shells were very massive — of order 10-20 M — in 



reservoir of kinetic energy ( 


Smith et al. 


2007, 


2008a. 


2010: 


Smith & McCrav 


2007; 


Wooslev. Blinnikov & Heger 


2007 


)• 



similar to that of 



Smith fc McCravl ( 2007 ) have argued ba sed on a simplified analytical model, 



Falk fc Arnettl ()1977l ). that the high luminosity and long 



duration of these SNe can be explained by a SN colliding with a very massive 
and initially opaque CSM shell. We explore this idea here in more detail 
with a variety of possible CSM environments using numerical simulations. We 
suggest that the presence and shape of circumstellar shells can be a powerful 
tool to constrain the evolution of the progenitors of Type Iln supernovae. We 
investigate how the mass, speed, and morphology of such shells can influence 
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the evolution of a SN lightcurve. We undertake a parameter study of SNe 
with different masses interacting with a selection of possible circumstellar 
shells, both spherical and bipolar. From these simulations we calculate thermal 
emission profiles and compare them in order to constrain how the physical 
properties of circumstellar nebulae can influence the SN lightcurve, and to 
constrain the efficiency of converting kinetic energy to light. 
Our calculations are simplified in the way we treat the cooling of and ra- 
diation from the shocked gas, which we approximate as optically thin ra- 
diative cooling; by necessity; our application of these results is therefore 
limited in scope. An important point to note is that our approach is to 
simulate a variety of hypothetical SNe to demonstrate trends in how the 
lightcurve responds to changes in SN and shell properties. We are not at- 
tempting a quantitative fit to the observed data for any individual SN. This 
has been purs ued for a few relatively nearby and well-observ e d SNe Iln, such 



as SN f988Z (Turatto et al 



19991) . SN 1994W 



1993 



Chugai et al 



Chugai fc Danziger 



1994; 



20041) . and SN 1998S 



Aretxaga et al 



Chugai 



2001 ). where 



the CSM properties were derived from fitting the observed light curves and 
spectra. Those authors inferred massive precursor shell ejections in the few 
years before core collapse, although the energy demands and required shell 
masses for these were not as extreme as for SNe 2006tf and 2006gy. Our work 
here builds upon these earlier studies. 

We explain our adopted initial conditions and the numerical method in §2 and 
§3, respectively, and in §4 we discuss some details of the shock interaction. 
In §5 we discuss how the resulting light curves depend on various parameters 
and in §6 we discuss shell velocities, and how these may help to interpret 
observations. Finally, in §7 we interpret our results in context with the most 
luminous SNe Iln, and in §8 we provide a summary. 

We include electronic datafiles containing the results of our simulations with 

this paper. The L dat files contain the total luminosity [erg/s] as a function 

of time [s]. The V dat files contain both the volume averaged and mass 

averaged velocity of the shocked gas [cm/s] as a function of time [s]. A small 
sample of these tables is provided in Appendices [A] and [Bj 



2 INITIAL CONDITIONS 



2.1 Supernova model 



In our si mulations, we begin with a co r e-collapse SN in free expa nsio n as de 



scribed by 



Chevalier fc Fransson 



|l992l ). lMatzner fc McKeel |l999l ) and 



Chevalier 
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Table 1. Simulation input parameters 



Name M sn E sn M shell V(6 = 0) A/ wind n t cnd dE/E sn vj inal [km/s] 

[ Mq] 10 51 ergs [M @ ] [km s" 1 ] [M yr" 1 ] Q [yr pre-SN] % 10 3 [km/s] 



O01 


30 


1 


N/A 


200 


10 -4 


0.0 


N/A 


0.05 


4.39 


O02 


30 


1 


N/A 


200 


IO -3 


0.0 


N/A 


0.3 


3.42 


O03 


30 


1 


N/A 


200 


io- 2 


0.0 


N/A 


1.5 


2.57 


O04 


30 


1 


N/A 


50 


IO" 4 


0.0 


N/A 


0.108 


3.75 


A00 


30 


1 


0.1 


200 


— 

10 -4 


0.0 


2 


0.8 


2.85 


A01 


30 


1 


1 


200 


IO" 4 


0.0 


2 


5.05 


2.34 


A02 


30 


1 


6 


200 


IO" 4 


0.0 


2 


18.7 


1.79 


A03 


30 


1 


10 


200 


IO" 4 


0.0 


2 


25.5 


1.59 


A04 


30 


1 


20 


200 


IO" 4 


0.0 


2 


36.5 


1.27 


A05 


30 


1 


10 


200 


IO" 3 


0.0 


2 


25.6 


1.51 


A06 


30 


1 


10 


200 


IO" 5 


0.0 


2 


25.3 


1.61 


A07 


30 


1 


10 


50 


IO" 3 


0.0 


2 


31.5 


1.30 


A08 


30 


1 


10 


500 


IO" 3 


0.0 


2 


16.9 


1.66 (at500 days) 


A09 


30 


1 


10 


50 


IO" 4 


0.0 


2 


31.6 


1.36 


A10 


30 


1 


10 


500 


1Q -4 


0.0 


2 


16.3 


1.80 (at500 days) 


All 


30 


1 


10 


50 


IO" 5 


0.0 


2 


31.6 


1.46 


A12 


30 


1 


10 


500 


IO" 5 


0.0 


2 


16.3 


1.80 (at 500 days) 


B01 


10 


1 


10 


50 


IO" 4 


0.0 


2 


54.5 


1.52 


B02 


10 


1 


10 


200 


io- 4 


0.0 


2 


48.2 


1.83 


B03 


10 


1 


10 


500 


IO" 4 


0.0 


2 


37.0 


2.12 (at 500 days) 


B04 


10 


1 


25 


200 


IO" 4 


0.0 


2 (At = 5yr) 


65.1 


1.08 (at 500 days) 


C01 


60 


1 


10 


50 


IO" 4 


0.0 


2 


19.7 


1.18 


C02 


60 


1 


10 


200 


IO" 4 


0.0 


2 


14.5 


1.30 


C03 


60 


1 


10 


500 


IO" 4 


0.0 


2 


7.55 


1.46 



D01 


10 


1 


10 


500 


IO" 4 


0.9 


2 


42.1 


N/A 


D02 


30 


1 


10 


500 


IO" 4 


0.9 


2 


20.4 


N/A 


D03 


60 


1 


10 


500 


IO" 4 


0.9 


2 


10.9 


N/A 


E01 


30 


1 


10 


500 


IO" 4 


0.0 


4 


14.7 


1.79 (at 500 days) 


E02 


30 


1 


10 


200 


IO" 4 


0.0 


10 


24.3 


1.62 (at 500 days) 


E03 


30 


1 


10 


500 


10~ 4 


0.0 


10 


13.9 


1.79 (at 1000 days) 


F01 


30 


0.5 


10 


200 


IO" 4 


0.0 


2 


22.3 


1.13 


F02 


30 


2 


10 


200 


IO" 4 


0.0 


2 


27.4 


2.19 


G01 


6 


1 


6 


200 


10~ 5 


0.0 


2 


56.4 


2.07 (at 100 days) 


HOI 


1 


1 


1 


200 


IO" 5 


0.0 


2 


42.7 


4.77 (at 100 days) 



(1) Measured at 250 days unless indicated otherwise. 



(2005), which gives a density profile divided in two segments: The inner part 
has p ~ r~ m , the outer part p ~ r~ b , with m = 1.06 and b = 11.7 for a 
progenitor star that still has a large hydrogen envelope at the moment of 
core collapse. The division between the two power laws lies at the transition 
velocity: 



t>t r = 3160* 



IE 51 



(5-m)(6-5) 
(3-m)(6-3) 

km s _1 



10 M 



(1) 



M„ 



(Chevalier fc Fransson 



1992; 



Chevalier 



20051 ). Using this profile we construct 
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three SNe, with different mass but equal energy. Because of the large value 
of b, the density drops very quickly at higher velocities. As a result, only a 
small fraction of the mass is moving fast, limiting the inertia. Therefore, the 
gas will slow down quickly when it collides with the circumstellar medium 
(CSM). Our standard massive-star SN has 30 M of ejecta mass and 10 51 erg 
of total kinetic energy, although we explored a range of SN ejection masses 
at 6, 10, 30 and 60 M Q , with total kinetic energies of 0.5, 1, and 2xl0 51 erg. 
We start each simulation of the CSM interaction at the moment were the 
supernova has expanded to 1 AU. Typical maximum velocity for the initial 
supernova is about 30 000 km/s. However, at this velocity the density is very 
low and the maximum velocity is quickly reduced to about 10 000 km/s by 
the collision with the surrounding medium, once the simulation begins. The 
distance it must travel to collide with the shell depends on the shell parameters 
(see 5I2T21. 

Our simulations do not include the effect of photo-ionization, nor do we take 
into account the effect of energy injection from radioactive decay. Our cal- 
culations simulate only the expected luminosity generated by the SN-CSM 
interaction shock front; our simulated light curves do not include emission 
from the expanding SN photosphere powered by diffusion of shock-deposited 
energy or from radioactive decay. These may affect the light curve at lower 
CSM-intcraction luminosities or very early times before the shock overtakes 
much of the CSM shell. Note also, that our supernova model is strictly spher- 
ical. Non-spherical SN ejecta outflows would greatly increase the parameter 
space and require a more complex calculation. 



2.2 Circumstellar shell model 

For the circumstellar shells, we take a variety of shell proper ties, but we 



focus on models reminiscent of the environment of rj Carinae (jSmith 



Smith et al 



2006 



20031 ). as such CSM properties have been proposed for some lu- 
minous SN Iln. Namely, we adopt a stellar wind with high mass-loss rate (10 -5 
to 10 -3 Mq yr _1 ) and moderate velocity (few hundred km s _1 ) for the steady 
wind phase before and after shell ejection, plus an expanding shell with ex- 
tremely high density that was ejected in an intermittent outburst reminiscent 
of giant LBV eruptions, occuring a few years before the SN. The mass-loss 
rate and velocity of the wind before and after the shell ejection are assumed 
to be identical. The supernova will therefore first encounter a (relatively) low 
density wind, then a short stretch of high density material and then once 
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again the low density wind after it escapes the shell. We explore a large pa- 
rameter space, covering a wide range of possible shell masses, velocities, wind 
mass-loss rates and ages. 



We also investigate the effect of a bipolar shell, such as might be e . 



Dwarkadas fc Owocki 



2002; 



Owocki 



cctcd 



2005 



by a rapidly rotating st ar (e.g. 

Smith fc Townsen I2OO7L The bipolar shape follows th e gravitational darken- 



i ng mo del for the wind of a rotating star as described by 
I2OO2J): 



Dwarkadas fc Owocki 



^ = l-tt 2 sin 2 (tf), (2) 
M(0) K ' W 



with Q = lu/uj c , u> — y / g/R the rotational angular velocity of the star and 
lo c the Kepplerian angluar velocity. Observations h ave shown th at the bipolar 
shell of 77 Carinae, for example, follows this shape (Smith 2006). The latitudi- 
nal angle 9 equals zero at the pole and 90° at the equator. Note that this set of 
equations only applies for radiatively driven winds. Should the star approach 



critical rotation during an eruption, mass could be focussed to the ec 



uator , 



2003). 



forming a flattened equatorial structure as well (jSmith fc Townsend 
This is not accounted for in these equations. 

The total range of parameters in our simulations is listed in Table [U In all 
cases we assume that the shell ejection lasted two years, though we explore the 
effect of different shell cross-sections by varying the velocity. Wind velocity 
and shell velocity are assumed to be the same, allowing us to use an analytical 
description, rather than a numerical model, for the shell morphology. 

The second to last column in Table [T] gives the efficiency of converting shock 
kinetic energy into radiated luminosity as found in our simulations, based 
on the input kinetic energy and the integrated luminosity in the light curve. 
This is the maximum efficiency corresponding to the bolometric luminosity 
output. The efficiency in converting shock kinetic energy to visual light must 
be comparable to or less than this value. 

Finally, Table [T] shows the the velocity of the SN remnant after it has collided 
with the shell. We measure this velocity at a fixed point in time, except where 
indicated otherwise; these exceptions are necessary due to the nature of the 
circumstellar medium, which may require a longer time interval before the 
supernova has broken through the shell. Also, we don't list a final velocity 
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for those SNe that interact with bipolar shells, since for thise simulations the 
velocity is angle-dependent. 



3 NUMERICAL METHOD 



We use the ZEUS 3D code ( Stone fc Norman 



1992; 



Clark 



1996) for our sim- 



ulations. The grid is spherical and two-dimensional, with 500 gridcells along 
the radial axis and 100 gridcells along the asimuthal axis, covering a 90 de- 
gree angle from pole to equator. We have seeded both the initial supernova 
and the circumstellar nebula with small scale density fluctuations (5% for the 
supernova ejecta and 1% for the circumstellar medium). This ensures that the 
supernova breaks up the circumstellar shell upon collision. 



3.1 Grid evolution 

In order to achieve a high resolution at the collision between supernova and 
circumstellar medium, the size of the radial gridcells decreases with the ra- 
dius. This gives us the highest resolution at the outer boundary. Since we 
need to maintain this high resolution at the collision fro nt, we use the mov- 
ing grid option that is part of the ZEUS 3D code (see 



Whalen et al 



2008 



van Veelen et al 



2009). At the start of the simulation the freely expanding 
supernova fills the entire grid, with the exception of the outer radial bound- 
ary, which is set to an inflow boundary condition with the parameters of the 
circumstellar medium that the supernova is running into. At the end of each 
timestep the code finds the highest radial velocity within 50 radial gridcells 
of the outer boundary. Using the velocity in this cell as a basis all gridcells 
are moved outward as well, with velocities: 

^ - ^rtrwr (4) 

with r[i] the radius of the gridcell with index i, which runs from to 500, ic 
the index of the gridcell in which the radial velocity is highest and v[ic] the 
highest radial velocity within 50 gridcells of the outer boundary. The physical 
conditions at the outer boundary are updated each time the grid expands to 
conform to the values of the circumstellar medium at that particular radius. 
In this way the entire grid is stretched in the radial direction, ensuring that 
a) the supernova remnant can never overrun the outer boundary; b) a high 
resolution is always maintained close to the outer boundary where the collision 
takes place and c) the inner boundary is fixed and does not move. (N.B. This 
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method works well as long as one has to deal with a strong shock. It is not 
recommended for subsonic expansion). 

A drawback of this method is that the circumstellar nebula is supposed to 
be static during the SN expansion, whereas speeds of the pre-shock CSM for 
luminous SNe Iln seen in narrow P Cygn i absorption f eatures tend to be o f 



order 100-500 km s" 1 (Smith et al 



2007 



2008a, 



2010: 



Trundle et al 



2008). 



(The CSM speeds listed in Table 1 essentially determine the radii of the shells 
and therefore their density for an assumed total mass.) However, the velocities 
of the SN ejecta expansion are much faster than those in the CSM nebula, 
such that any evolution of the nebula during the SN expansion phase can be 
considered small. The inner radial boundary is fixed at r = and does not 
move when the grid expands. The inner radial boundary and both azimuthal 
boundaries are set to reflecting boundary conditions so no matter can escape 
from the system. 



3.2 Radiative cooling 

In order to obtain a lightcurve from our simulation we in clude the effect of op 



MacDonald fc Bailey 



tically thin radiative cooling, using the cooling curve from 
(|198ll ). We extend this cooling curve to temperatures above 10 10 K by assum- 
ing that for these temperatures the cooling curve depends on the temperature 
as A(T) ~ \/T (i.e., Brehmstrahlung). 

Rather than use the cooling routine that comes as par t of the ZEUS 3D code, 



we implement a new numerical method, described by iTownsendl ([20091 ). This 
method uses exact integration of the radiative cooling function rather than 
the traditional implicit or explicit schemes. It is faster, more accurate and 
avoids the potential instability of the old radiative cooling method used in the 
ZEUS 3D code, which uses a Newton-Raphson implicit calculation scheme. 
The assumption of optically thin radiative cooling to generate our lightcurves 
has some drawbacks. The circumstellar shells used in our simulations have high 
densities and are therefore likely to be optically thick to Thomson scattering if 
fully ionized. However, at such high densities, it is difficult for the material to 
remain fully ionized because of fast recombination rates, we believe that our 
assumption is acceptable for our limited pruposes, at least as far as radiation 
in the optical part of the spectrum is concerned. The high density of the 
circumstellar shells makes it unlikely that ultraviolet the light from the SN 
itself can fully ionize them. Those areas of the shell that become photo- ionized 
will undergo recombination on a very short timescale. Typical mass density 
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in the shell is about 5 x 10~ 13 g/cm 3 (see figs. [T] through [S]). Assuming pure 
ionized hydrogen for the sake of simpl i city t his gives us an electron density n e 
of 6 x 10 11 cm~ 3 . lDvson fc Williams! (|1997l ) give of recombination rate of 



Nr = n%(T e ) (5) 

with /3 2 {T e ) = 2x 10 _10 T e _3/4 cm 3 /s. For an electron temperature T e of 
10 000 K, this gives us a recombination rate of 7.2 x 10 10 per second. So even 
if fully ionized initially, the shell will recombine very quickly compared to 
the expansion of thcshell. It takes the SN at least several days to reach the 
inner edge of the shell, so the effects of the initial ionization will most likely 
have disappeared by then, leaving only the remaining radiation from the ex- 
panding shock to photoionize the shell. This greatly reduces the number of 
free electrons that are available for scattering. Furthermore, although we use a 
shell with a smooth density structure (apart from the small random variations 
mentioned above), in reality circumstellar shells show a far more complicated 
structure of high density filaments interspaced between low density areas. 
Under these cir cumstances, the photons will tend to escape through the low 



density regions (jOwocki et al.ll2004t lOwocki fc Cohenll2007l) . Finally, scatter 



ing by itself does not necessarily change the shape of the emerging lightcurve 
since a photon can escape with little modification even after multiple scat- 
terings. Therefore, even though the electron scattering optical depth of our 
denser shells (under assmption of full ionization), can be as large as r e ^ 100, 
the true optical depth will be much smaller due to efficient recombination. 

The shape of the lightcurve will change if the diffusion time for photons to 
escape from the circumstellar s hell gets close to the actual expansion time of 
the SN ( Smith fc McCravLoOT ). However, this is only likely to affect the light 
curves at early times; the net effect would be a slower rise time to peak lumi- 
nosity and possibly a smoother peak, while diffusion is unlikely to substan- 
tially affect the overall efficiency of converting kinetic energy into radiation. 
The typical diffusion timescale of a photon through the shell is idi// = tD/c, 
with t the optical depth, D the cross-section of the shell and c the speed of 
light, whereas the expansion velocity is t cxp — D/V. With the expansion ve- 
locity V typically below 2000 km s _1 (see the shell velocity plotted in fig. [531 
and also typical final velocities in Table [T]) and lower for the denser, more 
optically thick shells, the photons have time to escape from the shell ahead of 
the expanding supernova. 

Although we have attempted to account for radiative cooling in a realistic way 



Supernova collisions with circumstellar shells 11 




R[AU] 

Figure 1. Temperature (left) and density (right) for simulation A03 at t=11.5 days after the start of the simulation. The 
supernova has not yet collided with the circumstellar shell. The front of the supernova expansion (~ 50 AU) is clearly visible 
because of the high (~ 10* K) local temperature. The small insets show how instabilities form in the thin supernova shell. 
Clearly, we are at the limit of what can be achieved with this grid-resolution. 

in our calculations, this is a difficult problem and our method is simplified and 
necessarily limited. Therefore, when interpreting our results, we concentrate 
on relative changes from one model to the next as we vary input parameters 
like mass and speed, rather than the absolute values of the luminosity for any 
individual model. As noted earlier, it is not our goal here to fit the observed 
lightcurve and derive corresponding physical parameters for any individual 
SN, but rather, we aim to understand how the variety of possible observed 
properties arises from different input parameters. 

4 SN-CSM INTERACTION 

Here we describe the general properties of our SN-CSM interaction simula- 
tions. As a SN interacts with the CSM, we observe three phases dictated by 
our assumed input geometry: (1) A fast interaction between the SN and the 
(relatively) low-density wind inside the shell, (2) a slower interaction as the SN 
shock pushes into the much denser medium of the massive shell, and (3) the 
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Figure 2. Similar to fig. [T] but at t=23 days. The supernova ejecta have reached the circumstellar shell. Note that the high 
temperature region has become extremely narrow. This is due to the high density of the shocked gas, which allows it to cool 
very rapidly. The temperature of the shocked gas increases, as more kinetic energy is converted to thermal energy. Again, the 
small figures show details of the supernova shell, which is extremely thin. Local instabilities are small. 



final expansion phase as the SN has broken through the shell and continues 
to interact with the wind outside the dense shell. 

To demonstrate the strong interaction between a core-collapse SN and a cir- 
cumstellar shell, consider simulation A03 (Table [lj , in which a 30 M© explo- 
sion collides with a 10 M circumstellar shell moving at 200 km s _1 . Figures [1] 
to S] show snapshots of the temperature and density of the expanding SN as 
it interacts with the CSM (movies of our simulations are provided in the elec- 
tronic edition). The high post-shock density encountered because of the very 
massive CSM shells we use causes the radiative cooling to be extremely ef- 
ficient in these simulations, sometimes reducing the internal energy of even 
the shockedgas to the point where the temperature reaches a minimum value 
of 1,000 K|^| (This lower limit is a matter of numerical convenience that we 
impose upon the calculation.) Since the temperature difference before and af- 



1 The fact that our calculations cool to a temperature as low as 1,000 K with a standard cooling prescription has far-reaching 
implications for understan ding dust formation in CSM-inte raction SNe. This is not the topic of our study here, but recent 

2009 ) have demonstrated that new 



observations of SN 2006jc (Sm ith et al 



2008bh and 2005ip (Smith et al 



2009 



Fox et al 



dust grains are seen to condense in the post-shock gas at the same time when strong X-rays and high ionization emission lines 
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Figure 3. Similar to fig. [l] to [3] and , but at t=86.8 days. The supernova is approx. halfway through the circumstellar shell. The 
temperature of the shocked region is much reduced (to ~ 10 5 K), because the high density of the shell reduced the expansion 
velocity. The small figures show the thin high temperature layers on each side of the shell. The shell is not perfectly spherical, 
but the instabilities are extremely small as they are compressed between the expanding supernova and the high density material 
of the shell. 

ter cooling can be quite large, we show the temperature of the gas before the 
radiative cooling has been taken into account, which is more representative 
for the wavelength of the emitted radiation. This is the same temperature 
that we use to obtain the lightcurves in §[5J adjusted for adiabatic expansion. 
At first, the SN ejecta expand quickly as the forward shock encounters the 
stellar wind, creating a layer of hot (several times ~ 10 8 K), shocked gas 
(fig. [I}. At the inner boundary of this high temperature zone (the reverse 
shock), SN material piles up and creates a shell. The interaction is (nearly) 
energy conserving at this point. When gas crosses the reverse shock, the kinetic 
energy of the SN is converted to internal energy and heats the shocked gas, 
pushing the forward shock into the CSM. The shell, which is very thin due 
to radiative cooling, is subject to thin-shell instabilities. However, these take 
time to form and the shell is expanding rapidly, which limits their opportunity 
to grow. As a result, the shell retains its basically spherical shape. Because 

are seen. With efficient cooling in the dense shock leading to the low temperatures in our simulations, dust formation may be 
a natural consequence. 
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Figure 4. Similar to figs.[T]through[3] but at t=173.6 days. The shock heated layer remains extremely thin, indicating a nearly 
isothermal shock. The shock temperature has decreases because sweeping up the circumstellar matter slows down the supernova 
expansion. The instabilities are somewhat larger now, but remain small compared to the overal scale of the expansion. 

only a small fraction of the SN material has a high velocity (see § 12. ip , the 
blast wave slows down quickly when it sweeps up the wind. This effect is 
greatly increased in the next phase when the SN ejecta collide with the dense 
circumstellar shell. 

Initially, the collision between SN ejecta and a massive circumstellar shell 
causes a rapid decrease of the forward shock velocity. This deceleration drains 
energy from the forward shock, and powers the main peak of the light curve. 
The density at the forward shock increases sharply as the shock overtakes 
more of the massive circumstellar shell. The layer of hot, shocked gas is com- 
pressed as the reverse shock starts to overtake the forward shock, which leads 
to an increase in local temperature (~ 10 9 K in fig. [2|. The high tempera- 
ture, combined with the high density of the gas makes the radiative cooling 
efficient. Therefore the thermal pressure of the shocked gas does not increase 
further. This, combined with the compression between the two shocks causes 
the hot gas layer to become quite thin and marks the transition from an energy 
conserving shock to a momentum conserving one. 

As the forward motion slows, the shock temperature decreases. The cool- 
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ing remains efficient, so the high temperature region, which is now at about 
< 10 5 K, remains thin (see fig. [3]). The thin shocked gas layer is subject to 
radiative cooling instability (the higher density regions cool more efficiently, 
leading to a loss of thermal pressure, which in turn leads them to be com- 
pressed to even higher density) . This can be observed in fig. [3] as variations in 
the local temperature in the shocked gas. However, like in the initial phase, 
the expansion of the SN occurs at higher velocity than the formation of the 
instabilities. Also, the shocked gas layer is compressed between two areas with 
very high density (the shell on the outside and the rest of the supernova on 
the inside), which inhibits expansion apart from the bulk motion of the shock. 
Therefore, there is no significant departure from spherical symmetry. 
Once the SN breaks through the shell, the forward shock may accelerate again 
due to the transition to much lower densities in the wind, though it will never 
reach the original high velocity because a large amount of energy has been lost 
to radiation during the shell collision phase. Also, the velocity of the unshocked 
SN ejecta piling up at the reverse shock decreases over time, limiting the shell's 
ability to accelerate in this later phase. As a result, the temperature of the 
hottest shocked gas is now limited to a few times 10 6 K (fig. |4j. The lower 
density at the interaction front makes the radiative cooling less efficient, which 
allows the hot gas layer to build up, though it never reaches its original size. 

4.1 Bipolar nebulae 

The collision between a SN and a bipolar nebula shows the same general 
pattern as described above, but is somewhat more complicated and modified 
by the shell geometry. Figures [5] through [5] show the same time frames for 
simulation D02, which models a collision between the same SN as in simulation 
A03 with a 10 M Q shell, but here the shell is bipolar. 

Initially (fig.E]) the simulations look the same as before, but they diverge once 
the SN hits the circumstellar shell. This occurs first at the pinched equatorial 
waist of the nebula, where the shell radius is smallest (fig. [6J . The collision 
squeezes the region of shocked gas into a very thin layer. In polar directions, 
the SN still expands into a lower-density wind and the hot gas layer remains 
wide. At a later stage, the interaction with the shell has slowed the expansion 
at the equator, leading to a lower shock temperature (~10 4 K in fig. [7]), 
whereas the shock temperature at the pole is still high because the shock has 
only just reached the circumstellar shell and has not swept-up enough mass 
to decelerate yet. 
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Figure 5. Temperature (left) and density (right) for simulation D02 at t=11.5 days after the start of the simulation. This is 
the equivalent of fig. [T] but with a bipolar nebula. At this point in time the supernova expansion is almost identical to the 
expansion in a spherical CSM. 
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Figure 6. Similar to fig. [2] (same timestep), but for simulation D02. At the equator the supernova has reached the shell and 
has been slowed down abruptly at the pole the supernova is still expanding into the wind. Note the difference in the hot gas 
layer, which has been squeezed by the collision. 
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Figure 7. Similar to fig. \3\ (same timestep), but for simulation D02. At the equator the supernova is moving through the shell. 
The slow shock has reduced the temperature of the hot gas zone to about 10 4 K. At the pole the supernova has finally reached 
the circumstellar shell. There the hot gas is still at a high temperature (10 7 K). 
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Figure 8. Similar to figs. [4] (same timestep), but for simulation D02. The supernova has broken through the shell at the 
equator, but is still plowing through it at the pole. The temperature at the equator is now high again (~ 10 K), whereas the 
shock over the poel has slowed down, lowering the local shocked gas temperature. 
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Eventually, the SN will start to break out of the shell at the equator first, 
while it is still inside the shell at the pole (see fig. [8J. When this happens, 
the shock at the equator will reheat to about 10 6 K, while the temperature 
at the pole remains low. Since the circumstellar shell has most of its mass 
concentrated at the pole (where it also has the largest solid angle), it takes 
much longer for the SN to break out in that direction. As a result, the shock in 
the polar direction will always be less energetic afterwards than at the equator. 
Because of the different times when the shock hits the equatorial and polar 
regions of the shell, different shock temperatures can be seen simultaneously 
We therefore suggest that simultaneous observations of multiwavelength (i.e. 
X-ray and visual) light curves may provide a way to distinguish bipolar from 
spherical shells, as we describe in more detail later. 

5 SUPERNOVA LIGHTCURVES 
5.1 General Properties 

The assumption of optically thin cooling, though a reasonable approxima- 
tion in optical wavelengths, breaks down for high frequencies. Most likely, for 
massive shell collisions, the early-time X-rays and UV would be completely 
self-absorbed and reprocessed into visual-wavelength luminosity. Therefore, 
rather than attempt to plot the emission as a function of the gas temper- 
ature, we concentrate on the bolometric luminosity light curves as a likely 
proxy for the visual lightcurves in later sections; this assumption may break 
down at late times when the shock becomes optically thin and X-rays can 
escape (see below). In addition, it is important to note that our light curves 
correspond only to radiative energy losses from the post-shock gas. We do not 
include the photospheric emission from the underlying SN itself, which could, 
in principle, be any type of SN. (It is the shell collision that leads to a Type 
Iln spectrum and the enhanced luminosity, rather than any intrinsic property 
of the SN.) 

As is shown below, the overall shape of the lightcurve for any SN-CSM col- 
lision model has the same general properties. Initially, the SN expands into 
the (relatively) low-density wind, starting at high luminosity due to its high 
velocity. As the expanding shock sweeps up more wind material, the expan- 
sion speed is reduced and the lightcurve shows a corresponding decrease in 
luminosity. 

Note that the behavior in this early phase — while the shock propagates 
through the wind on its way to reach the inner radius of the dense shell — 
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depends strongly on our assumed value for the inner radius of the shell and 
on the assumed time before the SN when the shell ejection finished {t en d in 
Table 1). If the SN had occurred immediately after the shell ejection stopped 
or while it was still in progress, then this early phase would not exist. This 
may be an important consideration in determining the early light curve shape: 
some luminous S Ne like SN 2006gy and SN 2005gj show a long and slow rise to 



peak luminosity (jSmith et al 



2007 



Prieto 



20071 ). while others are discovered 



at peak and decline immediately, as in the cases of SNe 1998S , 1997cy, and 



2006tf (Leonard et al 



2002 



Germany 



2000; 



Smith et al 



2008al) , suggesting a 



very rapid initial rise time. 

When the supernova reaches the circumstellar shell, which takes on the order 
of 10-25 days in most of our simulations, the expansion decelerates abruptly. 
This shows up as a rapid increase of the emission, because the fraction of 
kinetic energy converted into thermal energy is now high Also, the very high 
density of the shocked gas causes it to radiate very efficiently. As the SN plows 
through the shell, the emission decreases again due to the general decrease 
in shock velocity, but remains high compared to the emission from the initial 
phase. 

Once the SN has overtaken the massive shell and begins to expand into the 
outer low-density wind (at t ~140 days), the total emission decreases because 
the density of the gas that the SN collides with has decreased. Unlike the 
previous phases while the blast wave was expanding inside and through the 
massive shell, the densities are relatively low, and so optical depth effects are 
less likely to cause complete self-absorption of high energy photons. Thus, 
once the blast wave has broken through the outer boundary of a hypothet- 
ical massive shell, we would predict that soft X-ray emission could in fact 
be observed. Mass-loss rates derived from observations of this X-ray emis- 
sion would trace the normal wind mass-loss rate of the progenitor star in 
the years before it ejected the massive dense shell that led to the enhanced 
optical luminosity; meanwhile, the optical luminosity is still being emitted 
by teh dense shell. Therefore, one would not necessarily expect agreement in 
mass-loss rates de rived from observed optical and X-ray luminosities (see e.g., 



Smith et al 



20071 ). As the SN blast wave continues to expand into the wind, 



it gradually decelerates because the amount of swept-up gas increases over 
time. This leads to a steady and slow reduction in total emission in the years 
after the initial collision. 

Since both the circumstellar shell and the SN are spherically symmetric, the 
collision happens at the same moment everywhere. Similarly, the SN will break 
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through the shell at the same time all around its circumference. As a result, 
both sides of the main peak in the lightcurve have very steep slopes, and the 
change in X-ray emission would most likely be quite sudden. This is partly a 
result of our prescribed geometry of the shell, with a clean inner and outer 
boundary. Real circumstellar shells can show a wide variety of different ge- 
ometries, including multiple shells and high degrees of clumping, which can 
vastly change the appearance of the lightcurves. As one example, we explore 
the influence that a bipolar shape has on the emergent light curve. Our point 
here is not to provide an exhaustive grid of simulations of possible light curves, 
but to simply illustrate the behavior as we vary the parameters of the colli- 
sion in order to guide the interpretation of lightcurves of luminous SNe. The 
responses of the light curves to various parameters of the wind and shell are 
described in the following sections. 

The most important consequence of the SN-shell collision is that SN kinetic 
energy is converted to thermal energy and then lost to radiation. The efficiency 
of this conversion is a key parameter for interpreting the energy budgets of 
SNe Iln. For each simulation discussed below, we list the total efficiency in 
converting kinetic energy into radiated energy over the course of the simu- 
lation, E ra d/ Esn , in the second to final column of Table 1. We find a large 
range in the conversion efficiency, depending on the mass of the shell as well 
as the mass of the SN. For a circumstellar shell mass of 10 M , the efficiency 
is typically 15-30%. Efficiency increases with increasing density of the cir- 
cumstellar shell (higher shell mass, slower velocity, or both). The efficiency 
also increases for lighter SNe (higher ratios of M s h e ii/MsN), because of mo- 
mentum conservation and the greater deceleration of the fast SN ejecta. We 
elaborate on these points for individual cases below. 



5.2 No Shell, Just a Wind 



Since we are interested in investigating the effects of various properties of 
massive circumstellar shells, one might first ask what the collision looks like 
when there is no shell — i.e., when it is simply a collision between the SN 
and a dense steady wind. Figure |H] shows the bolometric luminosity emission 
lightcurve and the shock velocity (See §6) for a simulation where the circum- 
stellar medium contains no shell (O01), but just a dense wind with M = 10~ 4 
Mq yr _1 expanding at a speed of 200 km s^ 1 as one might expect for a 



massive luminous blue variable (LBV) progenitor ([Smith et al 



2007 



2010: 
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Figure 9. Lightcurve and reverse shock velocity for a SN expanding into a circumstellar medium that contains only wind 
(simulation O01). 



Trundle et al 



2008). Both the luminosity and post-shock shell velocity start 



high, but decline quickly as the SN sweeps into the dense wind. 
However, an important point to take away from simulation O01 is that the 
peak luminosity at early times is less than 10 41 ergs s _1 , and is therefore likely 
to be dwarfed by much stronger emission from an underlying SN photosphere 
(not shown in Fig. EH . A norm al Type II-P supernova, such as SN 1999em 
(Figure [TO] 



Leonard et al 



2002), has a luminosity during a ~110 day plateau 



of ~10 42 ergs s" 1 . This is lOOx stronger than the day 100 luminosity in simu- 
lation O01. Even SN 2005ip, which represents the lower end of luminosities for 
Type Iln core collapse SNe, h ad a late-time luminosity due to c ircumstellar 



interaction of 10 415 ergs s 1 ([Smith et al 



2009; 



Fox et al 



200J). The main 



consequence is that the more luminous class of SNe Iln require massive cir- 
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Figure 10. Lightcurves for four simulations without circumstcllar shells. (001-004) In most cases the luminosity fal ls well below 
peak values for a typical SN II-P photosphere (~10 42 ergs), illustrated by the observed lightcurve of SN 1999em j Leonard et al.l 



cumstellar shells, ejected in outbursts occuring shor tly before core collapse 



rather than steady wind s - as emphasized elsewhere (jSmith et al.ll2007 



Smith &: McCrav 



2008a; 



2007). This is also illustrated by a comparison between sim- 



ulation O01 and A00. these have the same input parameters except for a very 
low mass (0.1 Mq) circumstellar shell in the case of A00. Despite the low 
mass, the shell causes the total amount of energy converted to radiation to 
jump by more than an order of magnitude. 

Winds with higher density, either through high mass-loss rates (O02 and O03) 
or low velocity (O04) tend to produce higher luminosities through the colli- 
sion, as expected, but these enhancements are small compared to the effect 
of massive shells. (See also Table Q] for the percentage of energy converted 
into radiation.) The only "no-shell" simulation to produce a higher lumi- 



Supernova collisions with circumstellar shells 23 



nosity than that caused by even the smallest circumstellar shell is simula- 
tion O03, which assumes a mass-loss rate of 10~ 2 M yr _1 . Interestingly, this 
high wind mass-loss rate produces a late-time plateau with a l uminosity of 
~10 415 ergs s -1 , app ropriate for the late phases of SN 2005ip (jSmith et al 



2009; 



Fox et al 



20091) . A similar progenitor mass- loss rate was inferred for 



SN 2005gl, which had an LB V-like progenitor identified in pre-explosion data 



(Gal- Yam fe; Leonard 



SN 1988Z (|Aretxaea et al 



2009) .The shapes of these light curves also resemble 



1999), where the luminosity remained high for 



about a decade, indicating that the expanding supernova interacted with an 
extended circumstellar wind, rather than a sharply confined shell. Such a 
mass-loss rate is in excess of even the strongest L BV winds in their quies- 
cent states (i.e. 10 _3 M© yr _1 in the case of ij Car; 



Hillier et al 



(12001ft ). but 



i s comparable to smaller LBV eruptions like the 1600 AD event o f P Cygni 
( Smith fc Hartiganll2006 ) or the 1890 eruption of n Car ( Smith 2005 ). In other 



words, a steady "wind" with M = 10 _2 M© yr _1 is essentially the same as a 
sustained eruption (i.e. the total mass swept up by the shock is comparable). 
This is also the only simulation without a shell for which the radiative lumi- 
nosity exceeds values typically expected from the SN photosphere (Fig. [T0|) . 
One can expect that steady winds or sustained eruptive phases with even 
higher mass-loss rates or slower wind speeds will result in long-lasting light 
curves shaped like those in Figure I10| but with even higher luminosity. 



5.3 Shell masses 

The next group of simulations in Table Q] (A00 to A04) explore the effect 
that the circumstellar shell mass has on the evolution of the SN light curve. 
Fig. [TT] shows the total bolometric radiative luminosity for each simulation as 
a function time, compared to the light curve of a normal SN II-P. Because it 
takes more energy to break through a more massive shell, more kinetic energy 
is converted to thermal energy and then to radiative energy loss. Therefore, the 
higher the shell mass, the higher the luminosity peak. Also, it takes longer to 
break through a high mass shell, because the shock suffers more deceleration, 
so the duration of the peak luminosity will be longer for higher mass shells as 
well. It is noteworthy, that even the lowest mass shell (0.1 M©) simulation, 
A00, shows a clear peak and is therefore distinguishable from the pure wind 
interactions shown in fig. 1101 although in practice this peak might be lost amid 
the photospheric emission from the SN itself. 

All shells show sharp transitions at the beginning and the end of the main 
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Figure 11. Total bolometric luminosities for simulations AOO through A04 as a function of time. The higher shell masses 
cause higher luminosity peaks, for all other explosion and shell paramters held the same. Higher shell masses also cause more 
deceleration, so the shock takes longer to break through the shell, leading to a longer lasting peak in the lightcurve. Both the 
beginning and the end of the luminosity peak is marked by a sharp transition in all simulations, which results from our assumed 
inner and outer boundaries of the shell. The lightcurve of the SN II-P 1999em is shown again for comparison, as in Figure [lOl 



luminosity peak, but this is a direct result of our prescribed sharp inner and 
outer boundaries of the shells. It is a simplifying assumption and is motivated 
by the observed sharp outer boundary in s ome dense sh ells around massive 



stars, such as the Homunculus of r\ Carinae (|Smithll2006f ). but it is not neces- 
sarily true in all cases. It is likely that some objects will have smoother density 
transitions at the outer extent of the shell, and in those cases one expects the 
CSM-interaction luminosity to drop more gradually. The plateau is almost 
horizontal for the lower mass shells in our study, but changes to a shallow 
decrease with time for high mass shells. This decrease results from the fact 
that the high mass shells decelerate the blast wave to a greater extent as it 
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Figure 12. Total radiative emission lightcurves for circumstellar matter with different velocities (50 km s _1 in the left panel, 
200 km s _1 in the middle and 500 km s _1 on the right) and different wind mass loss rates (color coded lines). Obviously, the 
wind velocity has a major influence, since it determines how far the shell has travelled before the supernova hits it. Wind 
mass-loss rate makes very little difference, except in the very early stages. 



plows through the shell. The decrease in shock speed leads to a reduction in 
post-shock thermal energy and a lower emergent luminosity. 



5.4 Wind parameters 

In our simulations we vary both the wind velocity and mass-loss rate to ex- 
plore the influence of these parameters on the light curve. Figure [T^lshows the 
effect of the wind mass-loss rate and velocity on the bolometric SN lightcurve, 
by comparing simulations A07, A09 and All (left-hand panel in fig. [T2|) . which 
have identical parameters except for the wind mass-loss rate, which is 10 , 
1CP 4 and 10 -5 M Q yr _1 respectively. In these simulation the wind (and shell) 
velocity is fixed at 50 km s _1 . In the initial stage the difference is consider- 
able, as the higher density winds clearly create much stronger emission. Also, 
the high wind density in simulation A05 actually slows down the supernova 
expansion more than the other two, delaying the moment when the expansion 
reaches the shell, though not by a large amount. Since the shells are identical 
the lightcurves all have the same peak in the lightcurve. After the circumstel- 
lar shell has been swept up, the difference between the lightcurves is difficult 
to see. The 10 M shell slows down the SN expansion to such an extent that 
the effect of the wind mass-loss rate becomes negligible. Still, after more time 
passes the curves start to diverge, albeit slowly, with once again the highest 
mass-loss rate creating the highest emission. 

The middle and right-hand panels of fig. [12] show the same phenomena, but 
for wind (and shell) velocities of 200 and 500 km s _1 respectively. The results 
follow the same pattern. However, due to the higher velocties, the densities 




Figure 13. Total radiative emission lightcurves for collsions with different supernova mass and wind velocities. The higher the 
supernova mass, the lower the emission, as the high mass supernova has a relatively low velocity. 



are generally lower. As a result, the influence of the wind mass-loss rates in 
the final stages is lower for the simulations with wind velocity of 500 km s . 
Comparing the three panels of fig. [T2l shows the effect of wind velocity on the 
lightcurve. Obviously, lower wind velocities mean that the shell is closer to the 
star when the supernova hits is, which means that the entire time-frame of the 
interaction gets shortened. Also, the density in the shell is higher (p ~ l/ u )> 
whereas the cross-section of the shell is smaller, leading to a higher, narrower 
peak in the luminosity. 



5.5 Supernova masses 

Letting three different SNe interact with the same circumstellar shell pro- 
duces the light curves shown in Figll3[ which shows the bolometric lightcurves 
for three different supernova masses (10, 30 and 60 M Q : colorcoded lines), 
colliding with three different circumstellar shells (velocities at 50, 200 and 
500km s _1 : left, center and righ respectively). All three circumstellar shells 
have the same mass of 10 Mq. 

These lightcurves show two characteristic patterns: Because the kinetic energy 
in the SN is the same for all three simulations, the lower mass SNe have 
higher initial velocities. As a result, the peak in luminosity that results from 
the collision between the SN expansion and the circumstellar shells occurs 
earlier, and the peak luminosity is higher because of the greater energy per 
unit mass that is lost to radiation when the material is decelerated. The low 
mass SNe have less momentum (mA09 v \t)$ = m BoiVg 01 and mAog = 3msoi 
so rriBoi v BQi = 1TIA09VA09/ V%) , so they slow down and give up their kinetic 
energy more quickly during the collision. As a result, the 10 Mq supernova 
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produces a lightcurve where the flat plateau in the lightcurve peak is sharply 
angled, rather than horizontal as for the higher mass supernovae. The slope of 
this plateau may therefore provide a useful diagnostic to constrain the mass 
and momentum ratios of the underlying SN and CSM shell. This same pattern 
can be seen in all three figures. The essential result is that relatively lower- 
mass SNe (i.e. faster SNe) have higher efficiency in converting kinetic energy 
into radiation, while more massive SNe have more momentum and therefore 
lose less of their kinetic energy to radiation. This exactly the opposite of the 
effect of the shell mass, which produces a higher efficiency when the shell 
is more massive. Therefore, the highest efficiency will be achieved for those 
collisions wherein a relatively low mass SN collides with a relatively high-mass 
shell. (See also fig. [24]). 
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The influence of the wind and shell velocity is similar to that observed in 
Fig. [12J Higher expansion speeds stretch out the duration of the light curve 
and lower the peak luminosity, because the collision takes place later and over 
a longer time, and the shock plows through a lower-density shell for the same 
shell mass. 

Extreme cases of low SN mass can be seen in fig. [HI which shows the lightcurves 
for simulations G01 and HOI, where supernovae of 6M and 1M respec- 
tively collide with shells of equal mass. The resulting lightcurves show peaks 

-|45 



most luminous SNe observed to date ( 


Ofek et al 


2007 


Ouimbv et al. 


2007 


Miller et al. 


2009: 


Gezari et al. 


2009: 


Smith et al] 


2010) 


. However, due to the 



low masses of the circumstellar shells, the bright peak fades quickly as the 
shell is swept up within just a few weeks, which is faster than the observed 
examples. 



5.6 Supernova Energy 

Altering the total energy of the initial SN explosion also can change the appar- 
ent shape and luminosity observed in the light curve during its col l ision with a 



circumstellar shell. Indeed, in the case of SN 2006gv. lSmith et al.l (|2010l ) mea- 
sure a total energy (E ra d + kinetic) of at least 5xl0 51 erg. Figure H"5l shows 
the lightcurves resulting from simulations F01, A03, and F02, where three 
different SNe of the same mass but kinetic energy of 0.5, 1, and 2 x 10 51 
erg, respectively, all collide with the identical circumstellar shell of 10 M 
expanding at the same speed of 200 km s _1 . 

In the discussion above, we found that higher SN ejecta speeds and lower SN 
masses (a result of assuming that they all have the same explosion energy of 
10 51 erg) was a key factor contributing to a high peak luminosity. The key 
ingredient of higher ejecta speeds can also be achieved with more total energy 
in an explosion, so we explored this as well. As one might naturaly expect, 
more energetic SNe lead to higher peak luminosities because they give up more 
of their initial energy as their faster ejecta suffer a sharper deceleration during 
the collision. The light curve peak is also narrower (shorter in duration) for 
the more energetic and faster SNe because it takes less time to overrun the 
same shell. 

The net effect of altering the SN energy is similar to that of changing the SN 
mass (but keeping the same energy), mainly because of the strong influence 
of the SN ejecta speed (i.e. compare Fig 14 to the middle panel of Fig 13). 
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Figure 15. The influence of SN kinetic energy on the lightcurve. The higher the total energy, the higher the peak in the 
lightcurve and the shorter its duration, due to the increase in velocity. 



Comparing F01, A03, and F02 in Figure 15 and Table 1, we see that SNe 
with higher explosion energy had higher peak luminosities, but also more 
total radiated energy and higher efficiency in converting shock energy into 
radiation, due to their higher speeds as discussed above. As we will see below, 
a major difference between these three SNe of different initial kinetic energy 
is seen in their final blast-wave speed after the shock overruns and exits the 
circumstellar shell, providing a potentially useful observational diagnostic (see 
§6). 



5.7 Bipolar nebulae 

So far, all our light curves have resulted from the collision between a spher- 
ical SN and a spherical circumstellar nebula. In contrast, Fig. [16] shows the 
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Figure 16. Bolometric lightcurves for the collision between three different supernovae with m = 10 (B01), 30 (A09) and 
60 (C01) Mq and a bipolar circumstellar shell. The Peaks in the lightcurve show the highest luminosity for the lowest mass 
supernova just as in fig. 1131 The peaks are much roundar than for the collisions between supernovae and spherical shells. 



bolometric light curve produced by the collision between the three SNe of 
three different masses and a 10 Mq bipolar nebula. Fig. Ql)] is analogous to 
Fig. [T31 but with a range of speeds in a single shell as a result of its bipolar 
geometry instead of a range of speeds in three different spherical shells. As 
with the spherical nebulae, the lower mass SNe tend to produce higher peak 
luminosity in the bipolar case because of their higher SN ejecta speeds. How- 
ever, unlike the collisions between SNe and spherical nebulae, the luminosity 
peaks have smooth curves an d more gradual slo pes, somewhat reminiscent of 



the lightcurve of SN 2006gy (jSmith et al 



2007 ). In our simulations, at least, 



this smoothness results from the bipolar shape of the nebulae. Rather than an 
instantaneous collision between the SNe and a circumstellar shell, the inter- 
action starts gradually, with the collision beginning first at the equator and 
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than eventually spreading to the pole. An analogous transition happens when 
the SN breaks out of the shell. Again, this happens first at the equator and 
only much later in the polar region. As a result, shocked gas regions with rad- 
ically different temperatures and densities can exist simultaneously, as shown 
in Figs. [5] through [S] One might imagine that a smoother light curve may 
also result from a smoother transition in density at the outer boundary of the 
shell. 

A side effect of this situation would be that the possible onset of X-ray emis- 
sion would be more gradual and not coincide with the drop in total luminosity 
as the SN breaks out of the shell. The X-ray curves are expected to be strong 
when the supernova collides with a wind rather than a shell, both due to 
higher shock velocity and lower optical depth. If the shell is spherical this 
transition happens everywhere at the same time. In the case of the bipolar 
shell, the supernova breaks out at the equator long before it can break out 
at the pole. Therefore, part of the shock may already generate observable 
high energy photons, while another part is still plowing slowly through the 
shell and emitting at much lower temperature with all high energy emission 
being absorbed. Again, this effect may be relevant to the well-studied case of 
SN 2006gy, where the progenitor mass-loss rate inferred from the observed X- 
ray emission and Ha luminosity is in severe disagreement with the mass-loss 
rate needed to pow er the continuum luminosity in a CSM interaction scenario 
( Smith et al From fig.[5]through[8l one might understand this apparent 

contradiction if, for example, the X-rays are generated at the equator where 
the forward shock has already broken through the shell, whereas in the polar 
region the shock is still plowing through the dense massive shell and thereby 
powering the continuum luminosity. 



5.8 Time frames 

If the time interval between shell ejection and the SN changes, this too will 
influence the shape of the light curve. We investigate this effect with simula- 
tions E01, E02, and E03, with the resulting light curve shown in Fig. [T7] If 
there is a longer interval of time between the precursor shell ejection and core 
collapse, there are two effects. First, a given shell can travel further away from 
the star and will therefore have a lower density for the same shell mass. This 
will reduce the peak luminosity resulting from the shock interaction. The sec- 
ond effect of a larger time lag between shell ejection and core collapse is that 
it delays the onset of the strong CSM interaction phase. This can potentially 
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Figure 17. Bolomctric lightcurves for the collision a 30 Mq supernova and three 10 Mq shells, ejected at different times and 
velocities. The parameters of simulations E01 and E02 have been chosen in such a way that the inner boundary of the shell is 
at the same position, though moving with different velocities. As a result, the lightcurves are nearly identical. Simulation A03 
shows the lightcurve that results from collision with a shell that is much farther away from the star. As a result the peak in 
the luminosity is much shallower. 



lead to a second light-curve peak if one also considers the initial rise and fall 
of photospheric emission of the underlying SN that we do not include here. 
For simulations E01 and E02, the shell velocities and ejection times have been 
chosen so that the inner boundary of the circumstellar shell is at the same 
position for both shells, though they are moving at different speeds. As a 
result, the lightcurves arc quite similar in onset and duration despite their 
different speeds. Due to the difference in shell velocity, the shell in E02 is 
denser than in E01, resulting in a higher luminosity peak, which, however, 
quickly disappears as the shock slows down. Generally, the lower density of 
these shells (p <~ 1/r 2 ), results in lower luminosity peaks more than 1 yr after 
core collapse, with edges that are less steep. They do not show the round peaks 
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Figure 18. Total luminosity for simulations B02 and B04, demonstrating the effect of having a shell ejection that lasted over 
a longer period of time, making a thicker shell. These simulations are identical at first, except that the CSM shell in B04 has 
a larger outer radius at the same density, and thus has a higher total shell mass and remains at high luminosity for a longer 
time as the blast wave plows through this additional material. Since the cross section of the shell is larger for B04 the peak in 
luminosity lasts much longer. 



observed for bipolar shells (TTf5)) . This effect is seen most clearly in lightcurve of 
simulation E03, which is extremely slow in its evolution, remaining luminous 
for several years. 

Another parameter is the outer boundary of the massive shell, determined in 
our simulations by the duration of the shell ejection episode and the speed of 
the shell. In all simulations discussed so far, the duration of the shell ejection 
phase was At=2 yr, and we varied the outer radius of the shell by adjusting 
the speed of the shell. However, the duration of the shell ejection can vary 
as well . The 19th century eruption o f r\ Carinae, for example, lasted about 
20 yr (jDavidson fc Humphreys! Il997l ). although th e mass eject ion seems to 



have been concentrated in a shorter time interval (Smith 



20061 ). If the shell 
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Figure 19. Shocked gas velocities for the same simulations as in fig. 1111 The velocity drops very abruptly when the supernova 
hits the circumstellar shell, then rises again as the supernova plows through the gas, then rises once more as it breaks out and 
makes the transition from a radiative shock back to an adaiabatic one. The final velocity depends clearly on the mass of the 
shell. 



ejection occurs at the same mass-loss rate but lasts longer, then the shell will 
be thicker and have a larger total mass. The main effect of this is that the 
main peak of the light curve would last longer at a comparable luminosity. 
This is indeed the case, as we show in Fig 1181 which compares simulations 
B02 and B04. This is different from the case mentioned above referring to the 
speed of the shell. If a larger outer radius and longer duration to the light 
curve result from a faster shell speed, then the density is lower for the same 
mass and the resulting luminosity will be much lower (compare simulations 
B01, B02, and B03). 
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Figure 20. Shocked gas velocities for the same simulations as in fig. 1121 The wind velocity changes the entire shape of the 
lightcurve, as it changes the location of the shell, whereas the wind mass loss rate is only important in the initial stages. 
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Figure 21. Shocked gas velocities for the same simulations as in fig. 1131 The high mass supernovae, which have relatively low 
velocity slow down less than the low mass supernovae, but not enough to reverse the shock velocities. 



6 SUPERNOVA SHELL VELOCITY 

An important observational parameter for SNe Iln, in addition to their ra- 
diative luminosity and total radiated energy measured from light curves, is 
the expansion speeds measured from line widths in spectra. In most SNe, 
the ejecta expansion speeds are inferred from P Cygni absorption features 



lying photosphere can be seen ( 


Chugai & Danziger 


1994 


Turatto et al. 


1993; 


Salamancal2002 


Smith et alJ2009 


). Often, however, the underlying SN photo- 



sphere is masked by the bright and possibly opaque emission from the dense 
shell of post-shock gas that powers t he excess luminosity in SNe Iln (e.g., 
Chug ai et al Looi l [Smith et aU2008al) . Fortunately, the dense shell of shocked 
gas that piles up at the contact discontinuity in the SN-CSM collision can be 
seen in the intermediate-width wings of the narrow Ha e mission lines, for ex- 
ample, and typically has a speed of a few 10 3 km s _1 (jChugai fc Danziger 
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Figure 22. Reverse shock velocities for the same simulations as in fig. 1 151 illustrating the effect of SN energy on v(t). The more 
rapid deceleration of F02 occurs simply because the fastest ejecta reach the inner boundary of the dense CSM shell sooner. 
The higher energy SN has higher speeds at later phases during the collision. 
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In order to estimate how the presence of a circumstellar shell influences the 
velocity of the post-shock gas seen in Ha emission, we plot the velocity of 
shocked g 8jS ctS 0/ function of time: specifically, we plot the mass-averaged 
radial velocity of the gas between the forward shock (i?2) and the reverse 
shock (Rl) 



Jm Jo r2 s i n (9)P v rdrd9 
Im Jo r2 sin(d)pdrd6 



(6) 



We choose this method to quantify the shocked gas velocity because it gives a 
good result both in the adiabatic and radiative shock regime. Mass averaging 
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Figure 23. velocity curves for the reverse shock, shocked gas flow past the reverse shock and highest density shocked gas (the 
shell) for simulation A03. 



rather than volume averaging is more realistic since the luminosity is highly 
density dependent, so high density areas would dominate the emission. 
The behavior of the shocked gas velocity during the SN-shell collision gen- 
erally proceeds as follows: Initially, the velocity drops exponentially, because 
the blast wave declerates while sweeping up the surrounding wind. When the 
forward shock hits the circumstellar shell it practically halts and the reverse 
shock velocity drops abruptly as the gas between the two shocks is com- 
pressed. After the initial collision the shocked gas velocity increases again as 
the forward shock recovers. However, the velocity is now much lower since 
the interaction has become radiative, so much of the available energy has al- 
ready been lost. Also, the forward shock is now moving through a much denser 
medium. Although the forward shock accelerates again as it breaks out of the 
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shell and runs down a steep density gradient, it never recovers its original 
velocity, as is true for the reverse shock. The forward shock interaction does 
become nearly adiabatic again (see §E]), so the velocity remains higher than 
during the collision with the shell. By this time the shell has gained substan- 
tial mass through sweeping up the surrounding medium. Therefore its forward 
momentum is high, and the velocity remains nearly constant for a long period 
of time because the outer wind has insufficient mass to decelerate it. 
The abrupt loss of forward velocity in the reverse shock shell is, in principle, 
a robust characteristic of SN-shell interaction. Whether or not it is actually 
observable, however, is unclear. If the density of the pre-shock CSM is high, 
as it needs to be in the case of the more luminous SNe Iln, then one might 



expect the pre-sh ock gas to be very optically thick ([Smith fc McCray 



Smith et al 



2007 



20101) and the emitting surface may be well outside the shock. 



In that case, the observed Ha line profile would be dominated by the narrow 
component from photoionized pre-sho ck gas (typical l y a few 10 2 km s -1 ) and 



2001 



Dessart et al 



2008) out 



broad electron scattering wings (e.g., IChugail 
to a few 10 3 km s _1 . This is indeed thou ght to be the case for SN 2006gy, as 



discussed in detail bv lSmith et al.1 (|2010l ). 

The final shocked gas velocity, on the other hand, should be easily observable 
in all cases because of lower optical depths at larger radii and at late times, 
and may therefore provide an unambiguous constraint on the CSM mass and 
SN energy. In Figure [19l which shows the velocity for the same simulations as 
fig [TlJ we can see that the final velocity does in fact depend strongly on the 
mass of the circumstellar shell. If the shell mass is relatively high, the velocity 
decreases by a larger amount as momentum is conserved. 
In fig. [201 we show the shocked gas velocity for the same simulations as Fig- 
ure [12l This demonstrates the effect of the wind velocity and mass-loss rate 
on the shocked gas velocity. As all SNe have the same total energy in these 
simulations, the higher mass SNe start out with lower velocity. As can be 
seen, the mass-loss rate only matters in the initial stage, before the collision 
between the SN and the circumstellar shell. The wind velocity does make a 
significant difference as it determines the location of the shell relative to the 
star and therefore the timetable of the interactions, but it does not strongly 
influence the final speed of the shocked shell. 

In fig. [2lJ which shows the reverse shock velocity for the same simulations as 
FigureQSl we demonstrate the effect of the mass of the SN ejecta on the reverse 
shock velocity. As all SNe have the same total energy in these simulations, 
the higher mass SNe start out with lower velocity. As the ejecta collide with 
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the 10 M Q shell, the lower mass SNe slow down more, since they have less 
momentum. Still, after the collision, the lowest mass SN is still moving with 
the highest velocity and is the first to break out of the circumstellar shell. 
Finally, fig. 12121 shows the reverse shock velocities for SNe with the same mass, 
but different total energy. The high energy SNe start out with higher velocities, 
but also lose more energy in the collision. Since the collision takes place at 
an earlier stage, they slow down more, so the final difference in velocities is 
much less than initially. However, the total SN kinetic energy is perhaps the 
most influential factor in determining the final shell speed. 
Because of the specific nature of the collision, their are several different fea- 
tures that can determine the observed shock velocity of the SN. To illustrate 
this effect, we show three alternative velocity curves in Fig. [531 a U based on 
simulation A03. A) The velocity of the reverse shock, B) The velocity of the 
gas that has passed through the reverse shock, and C) The velocity of the 
shocked gas that has the highest density. Initially, all three curves move to- 
gether. The gas flow of shocked gas is slightly faster than the reverse shock 
itself, since this gas is actually moving through the shock. Before the collision 
with the circumstellar shell the highest density of shocked gas is at the re- 
verse shock, since this is where the SN ejecta piles up. This changes once the 
supernova collides with the shell. The reverse shock recoils from the collision, 
stopping completely or even reversing, depending on the density of the shell. 
The gas velocity inside the shock decreases as well, but not as much, since 
this is governed by the shock conditions. The shock changes from adiabatic 
to isothermal, restricting the velocity jump over the shock. As the SN then 
ploughs through the shell, the location of the highest density feature changes. 
It is no longer at the reverse shock, but rather at the forward shock, where 
gas from the shell is being swept up. Therefore, the flow speed of the high 
density feature actually becomes lower than the velocity of the reverse shock, 
since we are now sampling gas that is still in the process of being accelerated. 
Once the blast wave breaks out of the shell, the original situation is recreated, 
as once more the highest density occurs at the reverse shock and the shock 
conditions change back from isothermal to adiabatic. 

Examining Figures[l9]through[23l one can see that velocities measured in spec- 
tra obtained at early times can be powerful diagnostics of the rapid changes 
occurring during the initial shell collision, while later spectra that provide 
estimates of the final coasting velocity of the CDS are key diagnostics of the 
energy and momentum budget of the explosion. A potential complication for 
the early-time velocities, especially with more luminous SNe Iln, may arise of 



40 A. J. van Marie et al. 



the inner CSM is very optically thick. If the CSM outside the shock is highly 
opaque, then a radiative precurso r may cause the photosphere to reside out- 



side the shock (Smith et al 



2010I ). in which case the observed velocities are 



not indicative of the true expansion speed. In any case, combinations of pho- 
tometry and spectra at early times while the SN is still on the rise to peak 
are quite valuable in breaking the degeneracy of various models. 



7 DISCUSSION 

7.1 The Influence of SN and Shell Properties 

In the previous sections we have shown how circumstellar shells can influence 
the evolution of both the observed SN lightcurve and the observed velocity. 
The presence of a substantially dense circumstellar shell always causes an 
increase in the radiative luminosity, lasting until the blast wave breaks through 
the shell. The height of this luminosity peak depends primarily on the density 
of the shell (and so, also on its total mass and speed), in the sense that denser 
shells invariably lead to higher luminosities for the same underlying SNe. The 
duration of the luminosity peak is a direct consequence of the time it takes the 
blast wave to propagate though the shell, so it depends on the total mass of 
the shell, its expansion speed, and its inner and outer radii (i.e. the duration 
of the pre-SN ejection episode). A relatively more massive shell produces a 
slower blast wave, increasing the duration of the light curve peak and causing 
a higher luminosity. A faster expansion speed for the shell will also stretch the 
duration of the light curve peak by increasing its outer radius, but will make 
it less luminous for the same mass. 

In our simulations, typical luminosity peaks for spherical shells tend to have 
a flat plateau, which is either horizontal, or angled downward as the shock 
velocity decreases over time. The beginning and end of the light curve peaks 
are clearly defined with sharp edges, but this is just a result of our simplifying 
assumption that the shell has sharp inner and outer boundaries; real shells 
may have more complicated density profiles. These characteristics tend to 
disappear if the nebula is bipolar in shape, because different latitudes in the 
bipolar shell are hit by the blast wave at different times, and so the light curve 
shape is smoother. 

Whereas the total luminosity and the visual luminosity peak when the blast 
wave collides with the circumstellar shell, the temperature of the emitting 
gas decreases as the shock slows. A gradual decline in the characteristic tem- 
perature inferred from the continuum slope in visual-wavelength spectra or 
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multi-band photometry has been seen in several well-studied exa mples of very 



luminous SNe Iln, such a s SN 2005gj, SN 2006gy, and SN 2006tf dPrieto 



2007 



Smith et alJl2008aU2010l ). The X-ray drop may not be observed if initial phases 
are optically thick and X-rays are fully absorbed and reprocessed. Trapping 
at high optical depths is an effect that we have not included directly in our 
simulations; we consider it likely, therefore, that the visual radiation will trace 
the bolometric luminosity at early times, as we discussed earlier. This is why 
we have shown the bolometric luminosity light curve in our plots. Luminos- 
ity at high energies increases again once the SN breaks out of the shell and 
interacts with the (relatively) low-density wind outside the shell. This change 
in shock temperature is less clearly defined if the shell is bipolar, because 
both high- velocity and low-velocity interactions can occur simultaneously in 
different parts of the shell. A clumpy CSM may produce a similar effect. 



The observed velocity evolution of the dense post-shock H shell depends 
strongly on the CSM density and SN energy in our simulations. This velocity 
decreases steeply in the earliest phase of the expansion when the blast wave 
sweeps through the wind inside the dense CSM shell, and then it takes another 
drop when the shock hits the circumstellar shell. However, these velocities in 
the earliest phases may be difficult to observe because of high optical depth 
effects that are not taken into account in our simulations, as noted above 
for the early light curve shape. The characteristic velocity observed after the 
SN/shell collision ends depends on the shell mass, SN mass, and the total 
explosion energy, and is typically 1-3 x 10 3 km s _1 in our simulations. This 
i s comparable to th e observed line widths in lum inous SNe Iln like SN 2006tf 



( Smith et al 



km s" 1 in SN 2006gy (|Smith et al 



2008af ) o r SN 2005gi JPrietoll2007l) . The faster speeds of -4,000 



2010) imply a high er energy explosion and 



Smith et al 



a relatively high-mass SN. Indeed, 
energy of at least 5xl0 51 ergs for SN 2006gy. 



(2010) estimated an explosion 



The mass and initial speed of SN ejecta (and hence, the total explosion en- 
ergy) also influence the evolution of the velocity. SNe with higher ejecta mass 
have higher inertia and are decelerated less, but they also have slower initial 
expansion speeds for explosions assumed to have the same total kinetic energy, 
and so they can end up with slower final expansion speeds. A more energetic 
and relatively more massive SN explosion will emerge from the shell collision 
episode with a faster final shock speed. Since there is some degeneracy in any 
one type of observed property, spectral observations of the pre-shock CSM 
speed, the post-shock shell, and the SN ejecta speeds (if they can be seen) 
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Figure 24. The total efficiency (%) in converting shock kinetic energy into radiated energy from Table 1 for several represen- 
tative simulations, plotted as a function of the ratio of CSM shell mass to the total mass involved (SN ejecta + shell). The 
unfilled circles (and dotted line) represent our baseline simulations with SN ejecta with 10 51 erg running into 200 km/s shells 
of various masses ("A" models, plus B02 and B04). The unfilled triangles and squares are similar but for CSM speeds of 50 
and 500 km s , respectively. The filled circles are models F01, A03, and F02, showing the effect of different explosion energy 
for the same shell paramters. The X's show models G01, H01, and C01 (all with V eX p = 200 km s _1 ), special cases that have 
extended CSM mass or lower-mass SN ejecta. 



are valuable to combine with estimates of the luminosity from photometry to 
derive the physical properties of the CSM interaction. 



7.2 The Shock Conversion Efficiency 

Since some very luminous SNe Iln have measured values for their total radi- 
ated energy approaching or even exceeding the canonical SN explosion kinetic 
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energy of 10 51 erg, the efficiency at which they convert some fraction of their 
initial kinetic energy into post-s hock thermal energy a nd then radiation is 
key. In the shell-shocked model (jSmith fc McCravll2007l ). high efficienies are 
allowable because of the large radius at which shock energy is thermalized, 
allowing the SN to radiate before it expands and loses that thermal energy 
adiabatically. The second to last column of Table 1 lists the efficiency of this 
conversion as the ratio of the total energy lost via radiation in each simulation 
to the initial explosion kinetic energy of the SN, or dE/EgN- In Figure 12^1 we 
plot this efficiency as a function of another ratio, which is the CSM shell mass 
compared to the total mass in both the SN ejecta and CSM. We show the 
results for several simulations to demontrate various trends. 

The basic result is that efficient conversion of SN kinetic energy into radiation 
via CSM interaction requires a CSM mass that is comparable to or larger 
than the mass of the SN ejecta. The primary criterion for luminous SNe Iln 
that result from core-collapse SNe is therefore the presence of several M Q of 
circumstellar gas which must have been ejected very shortly before the SN. 
Explosions of very massive stars can have CSM interaction that is not very 
luminous if the CSM mass is small compared to the SN ejecta mass (as long 
as the SN ejecta are slow and heavy for a standard energy). A very effective 
way to convert a larger fraction of the total initial energy (more than half) 
into radiation is to have a more extended CSM shell at the same density, as in 
simulation B04, tracing mass loss for a longer time prior to the SN explosion. 
Of course, the longer a simulation runs into CSM material, the more kinetic 
energy can be converted into light — if one waits for ^100 yr or more, an 
extended SN remnant can tap a significant fraction of the total energy. Our 
aim here, however, is to study objects that do this very quickly in ~1 yr and 
thereby produce high luminosities during the initial light curve peak. 

The pre-shock CSM speed also has some minor effect on the efficiency, in the 
sense that slower CSM speeds lead to denser environments that trap more 
of the available kinetic energy because of their denser post-shock gas, and 
consequently, more efficient cooling. Also mildly influential is the speed of the 
SN ejecta, or equivalently, the SN explosion kinetic energy. More energetic 
explosions are more efficient in converting their available energy reservoir to 
radiation due to the higher velocity drop at the reverse shock. Thus, mild 
increases in explosion energy offer an alternative to exceedingly massive CSM 
shells in order to produce very luminous SNe Iln. Again, however, the CSM 
must be extended and massive in order to maintain that high luminosity for 
an extended time. 
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Time (Days) 

Figure 25 . This plot compares a few selected models to light curves of the two mo st luminous SNe Hn known. SN 2006gy 
(data from lSmith et al.l l|2007l . 120083 )1 is shown with unfilled circles, and SN 2006tf l|Smith et al.ll2008al ) is shows with filled 
circles. The models shows are A04 (30 Mq SN and 20 Mq shell; dotted grey), B02 (10 Mq SN and 10 Mq shell; solid grey), 
and B04 (10 Mq SN and 25 Mq shell; dark grey). Model B04 is the same as B02 except that the shell has a larger outer 
radius (it is thicker at the same density) and therefore has a larger total mass. While these models may account for the peak 
luminosities of SN 2006tf, they fall short of the peak luminosity for SN 2006gy and they fade too quickly for both. It is likely 
that simulations with more extended and more massive CSM shells or more energetic SNe need to be explored in these two 
particuar ca s es. Th e very late time data at around 400 d for SN 2006gy may have some contribution from a light echo (e.g., 
ISmith et all l|2008ch l. but the late-time luminosity after 1 yr for SN 2006tf is d ominated by strong ongoing CSM interaction 
because strong Ha emission is seen in the late-time spectrum l|Smith et al.ll2008alV For a typical Type Hp lightcurve, see Figs. 
[TOlandfTTI 

7.3 Comparisons to observations of luminous SNe Iln 

The central motivation for this study was to explore shock interaction with 
dense pre-SN CSM shells as a possible engine for the visual light from the 
emerging class of extremely luminous SNe, and to ask whether observed light 
curves can be compared to expectations of hydrodynamic simulations in or- 
der to constrain the physical properties of those shells and the underlying 
SNe. Below we briefly mention two recent well-observed examples that have 
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been our primar y objects for comparison: (1) SN 2006gy ([Smith et al 



Ofek et al 



2007 



20071 ) was the first of these super-luminous SNe that raised many 



questions about our understanding of the power so urces for these obje cts and 
about massive star evolution, and (2) SN 2006tf ([Smith et al.ll2008at ) which 
was nearly as luminous. Both have optical spectra of Type Iln suggesting 
the presence of dense CSM, although SN 2006tf appears to fit the canoni- 
cal picture of CSM interaction as the power source with fewer complications 



(Smith et al 



2008al) . Both appear to have high optical depths at early phases. 



The energy sources for these two SNe are of particular interest because, un- 
like more common SNe Iln at lower luminosity, their total radiated energy 
severely taxes the total SN energy budget: the energy radiated in visual light 
was 0.7xl0 51 ergs for SN 200 6tf ([Smith et alj|2008ah and 1.6xl0 51 ergs for 



SN 2006gy (jSmith et al 



120071) . Making a bolometric correction based on the 
observed temperature yields E rac i ~ 2.5 x 10 51 erg for SN 2006gy, and includ- 
ing the kinetic energy remaining i n the shell pushes the total initial explo- 



sion energy to at least 5x10 erg (jSmith et al 



2010(). For SN 2006gy, a pair 



instability SN or diffusion from an opaque shocked shell h ave also been sug- 



gested as possibilities for powering the observed luminosity ( Smith & McCrav 



2007; 


Smith et al. 


2007: 


Ofek et al. 


2007: 



Wooslev. Blinnikov fe Heger 



2007), 



although det ailed analysis of it s spectral evolution favors the opaque shocked 



shell model ( Smith et al 



2010). 



Two other SNe, SN 2005 ap (jOuimbv et al 



2009: 



Gezari et al 



20071) and SN 2008es ( Miller et al 



20091 ). have also been discovered recently to be among 



the most luminous SNe known. In fact, their peak luminosities are some- 
what higher than SN 2006gy, although they faded more quickly. We do not 
consider these for direct comparison with the same type of model discussed 
here because their spectra are not of Type Iln, but rather, they had normal 
broad-lines in their spectra indicating a photosphere receding through fast 
ejecta. The lack of narrow emission lines makes it likely that their radiation 
is produc ed primrily by diff u sion f rom an opaque shocked envelope, as in the 



model of 



Smith fc McCravl (|2007h . but with a smaller envelope mass than 



for SN 2006gy. The parameters in some of our simulations with slower (and 
therefore more dense) CSM and higher conversion efficiencies, such as B01, 
might be appropriate for these objects if diffusion were properly accounted 
for. Alternatively, it has re cently been suggested that these SNe m ay powered 



by the birth of magnetars (jKasen fc Bildsten 



2009; 



Wooslev 



20091) . energizing 



the opaque SN ejecta from within. 

We did not tune models specifically to fit the observed light curves of SNe 2006gy 
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and 2006tf, but we did explore a range of parameters for combinations of SNc 
and CSM shells comparable to relevant parameters estimated from observa- 



tions (Ofeketal 



2007 



Smith et al 



2007 



2008a . 



2010: 



Wooslev. Blinnikov fc Heger 



20071 ). Among our models, some of the highest peak luminosities were attained 
with relatively low mass (and therefore fast) SNe running into slower CSM 
shells, such as B01 and B02 (10 M SNe), or A07, A09, and All (slower 50 
km s~ 1 shells). Although these models achieved very high peak luminosities 
comparable to those of the most luminous observed SNe Iln, they faded too 
quickly, and so they fall far short of achieving the duration and total radiated 
light output of events like SNe 2006gy and 2006tf. Lower mass-SNe run out 
of momentum too quickly, or over-run the compact CSM shells too quickly. 
These low- and moderate-mass models may be applicable to SNe 2005ap and 
2008es, which attained high peak luminosities and faded quickly as mentioned 
above. Again, an important consideration is that our models do not include 
the possible delayed effects of diffusion when high optical depths are impor- 
tant, as one might expect for massive and slow (and therefore dense) CSM 



shells (Smith & McC 



rav 



curve (Falk fc Arnett 



20071 ). Including this may produce a smoother light 



19771 ). especially in the sense that it would dmpen and 



round-out the sharp initial peak in many of our simulations. Thus, it is possible 
that models such as B01 could be dominated by diffusion and may not appear 
as SNe of Type Iln, but confirming this conjecture requires additional work 
beyond the scope of this paper. Diffusion through an opaque shell, however, 
would not alter the later phases of our light curves after maximum light. 
There were some models that came close to matching the lightcurve behavior 
of SNe 2006gy and 2006tf with both high luminosity and relatively long dura- 
tions. These were models with SNe that had very massive and extended shells 
of 10-25 Mq, although even these seemed somewhat insufficie nt. Figur e 1251 
compares th e observed light curves of these two SNe (data from 



Smith et al.l 



2007 



2008al) to models A04, B02, and B04. All of these have shell expansion 



speeds of 200 km s 1 , c lose to the observed valu es of pre-shock material for 



SNe 2006tf and 2006gy (jSmith et al 



2007 



2008al ). 



7.3.1 SN 2006tf 

The light curve of SN 2006tf shows a slow and steady decline from peak 
luminosity, the approximate rate of which is reproduced in all three models 
shown. In our simulations, this decline rate is mainly due to the deceleration of 
the post-shock gas as the SN ejecta sacrifice energy to radiation. An interesting 
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result is that this decline rate during the main light curve peak for model B04 
roughly matches the radioactive decay rate of 56 Co, even though there is no 
luminosity from radioactivity included in these simulations — in other words, 
luminosity from shock-CSM interaction alone can in some cases mimic the 
radioactive decay rate. (One might also expect a similar decline from a steeper 
density gradient in the CSM shell, or a different velocity/density law in the 
SN ejecta. In SN 2006gy, for example, recent e vidence points to a Hubble-like 



expansion law in the CSM (jSmith et al 



2010). However, none of the models 



sustain the high luminosity for a long-enough time. SN 2006tf shows relatively 
high luminosity above 10 42 ergs s _1 even at very late times more than 1 yr 
after peak, consistent with a continuation of the same decay rate, but all three 
models in Figure [25] drop long before that time. In our simulations, this drop 
occurs when the forward shock exits the outer boundary of the dense shell 
and continues into the lower-density exterior wind shed by the star before it 
ejected the CSM shell. A similar drop in luminosity was observed in the light 
curve of SN 1994W, and w as also attributed t o the shock overruning the outer 
boundary of a CSM shell ([Chugai et al.ll2004l ). This sharp drop is not usually 
seen in the light curves of SNe Iln, however, suggesting that most SNe Iln 
have more extended CSM shells. 



We explored the effect that changing the outer shell boundary has on the light 
curves. Models B02 and B04, both shown in Figure l25l are identical up to the 
point when the radius of the outer shell boundary is reached in model B02. At 
this time, occurring around day 100, the luminosity in model B02 plummets 
as the shock runs out of dense CSM to interact with. In model B04, however, 
we simply continued the same shell properties to a larger radius by having the 
shell ejection occur with the same mass-loss rate over a longer time interval 
(At=5 yr instead of 2 yr, both ending 2 yr before core collapse). Thus, the 
slow decline from peak luminosity continued at roughly the same rate until 
day ^260, when its forward shock reached the outer shell boundary and the 
luminosity finally plummeted. This larger outer radius required a much larger 
shell mass, increased from 10 Mq in model B02 to 25 Mq in model B04. The 
general shape and luminosity of model B04 is similar to the 25 M q shell model 
that W oosley et al. (2007) suggested for SN 2006gy, although Ismith et al. 
(|2008ah noted that it also fit the early light curve of SN 2006tf well. A shell 
mass of 25 Mq is near the limit of what one might believe from a non-terminal 
stellar outburst if giant eruptions of LBV s like rj Carinae are representative 



(Smith et al 



2003; 



Smith fc Ferland 



20071 ). Shell masses beyond 25 M© also 
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begin to press the most basic limitat ions of even a very ma ssive star's mass 
budget at the end of its life (see, e.g. 



Smith fc Owocki 



2006) 



However, even the extremely massive shell of 25 M Q in model B04 cannot 
sustain a high luminosity long enough to account for the +1 yr observations of 
SN 2006tf (Fig. l25| because it drops too soon. Instead, the late-time luminosity 
of SN 2006tf seems to continue the same slow decline rate. The corresponding 
CSM shell mass that this would imply (roughly 50 M Q ) is staggering if the 
continued high luminosity were the result of simply extending the same shell 
to larger radii. One way to avoid such implausibly high shell mass would be to 
lower the density of the envelope but increase the total SN explosion energy 
above 10 51 ergs. Higher explosion energy leads to faster SN ejecta speed, 
so consequently, higher instantaneous luminosity can be achieved with lower 
shell densities. A larger explosion energy also relieves some of the strain on 
the efficiency of converting kinetic energy into light, since the total radiated 
energy of SN 2006tf is almost 10 51 ergs. 



7.3.2 SN 2006gy 



This SN presents addition al challenges , since the t otal radiated energy ac- 



tually exceeded 10 ergs ([Smith et al 



2007 



20101 ) . requiring a more ener- 



getic SN explosion no matter what the CSM properties are. None of our 
models were able to achieve the combination of the high peak luminosity 
and long duration of SN 2006gy, although our most energetic SN explo- 
sion was only 2xl0 51 ergs. Following the arguments above for SN 2006tf, 
then, one might expect that models with a more energetic explosion could 
matc h the light curve of SN 2006gy wi thout having implausibly high CSM 
mass (jWooslev. Blinnikov fc Hegerll2007T) . For example, the smooth light curve 
shape and slow rise to maximum are traits that were seen in our simulations 
with bipolar CSM shells, so one can imagine that a set of parameters similar to 
model D02 but with higher explosion energy may account for the light curve 
of SN 2006gy. This will be explored in a future paper. Diffusion from opaque 
shocked shell may also le ad to a smooth light curve appropriate for SN 2006gy 



(ISmith fc McCrav 



2007), and we have not included these high optical depth 



effects in our simulations. 



8 CONCLUSION 

In this paper we describe the influence of massive circumstellar shells on core- 
collapse SN lightcurves, with the primary motivation of trying to understand 
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the power source of extremely luminous Type Iln events and their relation 
to the diverse population of SNe Iln. We show how these circumstellar shells 
can indeed create extreme peaks in the luminosity such as have been observed 
in Type Iln supernovae like SN2006gy and SN2006tf. The luminosity of these 
SNe would require extreme amounts of 56 Ni if they are powered by radioactive 
decay, but if interactions in the CSM provide the power instead, then the 
shell masses and speeds that are required have reasonable precedent from 
obs erved properties of spat ially resolved shells around nearby massive stars 
(see 



Smith fc Owocki 



20061 and references therein) . 



Our investigation is by no means exhaustive. Pre-SN circumstellar shells may 
have a wide range of masses, expansion speeds, and radii, whereas we have 
adopted simplified shell geometries for illustrative examples. Additionally, the 
underlying SNe ejecta may have wide diversity in explosion energy, mass, and 
ejecta speed. In this preliminary investigation, our approach has been to vary 
each of these parameters individually to illustrate their influence on the light 
curve rather than attempting to accurately model any individual SN. We have 
attempted to find general ways to distinguish between different kinds of shells, 
using trends in the observed shapes of the lightcurves, their characteristic 
emission temperature, and observed shock speeds. We find that observations 
of the evolution of the shock speed is necessary to help break the degeneracy 



CSM help considerably as well (see e . g . , 



Salamanca 


2002; 


Smith et al. 


2007 



2008a 


2009 


2010 


Trundle et al. 


2008) 



20081 ) . This can be used to analyze the mass- 



loss history of massive stars in the last years prior to the explosion, which can 
be a powerful tool for studying the final stages of stellar evolution. Ultimately, 
we wish to know the physical origin of these SN-precursor events. 



The key result is that we confirm the large masses of ci rcumstellar shells 



hypothesized to account for some recent luminous SNe Iln ( 


Smith & McCrav 


2007; 


Smith et al. 


2007, 


2008a. 


2010; 


Wooslev. Blinnikov & Heeer 


2007 


), as 



well as the high mass and explosion energy of the underlying SNe. One can 
also produce a very high peak luminosity with lower mass if the shell is slow 
and the SN ejecta are fast, but a lower mass shell cannot yield both a high peak 
luminosity and a long duration of >100 days seen in some luminous SNe Iln. 
In fact, we suspect that even larger shell masses or larger explosion energies 
are needed to account for the observed light curves of the most luminous 
SNe Iln. Thus, more detailed attempts to model individual objects will be 
the focus of a second paper in this scries. 
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8.1 Future developements 



Further research is required for quantitative analyses of observed SN Iln 
lightcurves and to extract reliable absolute values of shell masses and SN ex- 
plosion energies. This must include adding the luminosity contribution from 
the underlying SN photosphere (powered by diffusion or radioactive decay) in 
cases where the CSM interaction luminosity is not extremely high compared 
to the ejecta photosphere, as well as using an improved treatment of post- 
shock cooling and radiative transfer at high optical depths in order to more 
accurately model the eme rgent radiation from the post- shock shells i n these 



simul ations. As noted by ISmith fc McCravl (|2007l ) and 



Smith et al 



(2008a, 



20101 ). it is likely that the CSM will be highly opaque, especially at the ear- 



liest phases, so the effects of radiative diffusion should be taken into account 
to properly model the emergent luminosity. Finally, all our simulations have 
adopted a Type II-P core-collapse SN density profile, but other types of SNc 
with different density profiles need to be investigated in a similar manner, 
since any type of SN can, in principle, be a Type Iln event if it runs into a 
dense H-rich environment. 
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APPENDIX A: LUMINOSITY 

This appendix contains a sample of our luminosity tables. The full tables can 
be found online. 
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Table Al. luminosity values for simulation O01 



time [sec] 


loglO(L) [erg/s] 


1.0000000E+05 


4.0854548E+01 


2.0000000E+05 


4.0681275E+01 


3.0000000E+05 


4.0674991E+01 


4.0000000E+05 


4.0714879E+01 


5.0000000E+05 


4.0557764E+01 


6.0000000E+05 


4.0541291E+01 


7.0000000E+05 


4.0630674E+01 


8.0000000E+05 


4.0439978E+01 


9.0000000E+05 


4.0538254E+01 


1.0000000E+06 


4.0346281E+01 



Table Bl. Shock velocity values for simulation O01 



time [sec] 


V(volume averaged) [cm/s] 


V(mass averaged) [cm/s] 


2.0000000E+05 


6.3296810E+08 


7.0761081E+08 


3.0000000E+05 


6.1563164E+08 


6.8375272E+08 


4.0000000E+05 


5.9150803E+08 


6.5441661E+08 


5.0000000E+05 


5.7608001E+08 


6.4455979E+08 


6.0000000E+05 


5.7573592E+08 


6.3470001E+08 


7.0000000E+05 


5.6127737E+08 


6.1689406E+08 


8.0000000E+05 


5.4696696E+08 


6.0784766E+08 


9.0000000E+05 


5.3620586E+08 


6.0075791E+08 


1.0000000E+06 


5.3759244E+08 


5.9993685E+08 



APPENDIX B: SHOCK VELOCITY 

This appendix contains a sample of our shock velocity tables. The full tables 
can be found online. 



